February 1, 2008 



LBL-38549 



^ ; Parton Model Sum RulesQ 

u 

P-I' 

Ian Hinchliffe and Axel Kwiatkowski 
Theoretical Physics Group 
O Ernest Orlando Lawrence Berkeley National Laboratory 

University of California 
■ Berkeley, California 94720 

O ' 

On 

Q-i, Abstract 

This review article discusses the experimental and theoretical sta- 
,JZ< • tus of various Parton Model sum rules. The basis of the sum rules in 

perturbative QCD is discussed. Their use in extracting the value of 
the strong coupling constant is evaluated and the failure of the naive 
s_! ■ version of some of these rules is assessed. 



1 This work was supported by the Director, Office of Energy Research, Office of High 
Energy and Nuclear Physics, Division of High Energy Physics of the U.S. Department of 
Energy under Contract DE-AC03-76SF00098. 



Disclaimer 



This document was prepared as an account of work sponsored by the United States 
Government. While this document is believed to contain correct information, neither the 
United States Government nor any agency thereof, nor The Regents of the University 
of California, nor any of their employees, makes any warranty, express or implied, or 
assumes any legal liability or responsibility for the accuracy, completeness, or usefulness 
of any information, apparatus, product, or process disclosed, or represents that its use 
would not infringe privately owned rights. Reference herein to any specific commercial 
products process, or service by its trade name, trademark, manufacturer, or otherwise, 
does not necessarily constitute or imply its endorsement, recommendation, or favoring by 
the United States Government or any agency thereof, or The Regents of the University of 
California. The views and opinions of authors expressed herein do not necessarily state or 
reflect those of the United States Government or any agency thereof, or The Regents of 
the University of California. 



Lawrence Berkeley Laboratory is an equal opportunity employer. 



ii 



1 Introduction 



One of the best tools to use in attempting to disentangle the structure of 
the nucleon is lepton-nucleon scattering where the lepton, whose couplings 
to electroweak gauge bosons is fully known, is used as a probe on the con- 
stituents of the nucleon. Lepton-nucleon scattering with large momentum 
transfer between the lepton and the nucleon is described in terms of the Par- 
ton Model [H, In its naive form this model describes the nucleon as a 
collection of non-interacting quarks and gluons. Lepton-nucleon scattering 
is then viewed as the sum of incoherent scatterings by the lepton off these 
partonic constituents. The description of these constituents is most con- 
veniently given in a frame where the nucleon has large momentum. If the 
nucleon mass is neglected, its momentum can be written as P^ = (p,p, 0, 0). 
A parton momentum can be written as Pf 1 = (zip, Zip,p t , 0) where pt ~ 300 
MeV is related to the scale of nucleon binding. A distribution function fi{z) 
is defined so that the probability that a parton of type % (for example an 
up quark) has momentum in the range P(z) to P(z + dz) is fi(z). The 
lepton-nucleon scattering rates are then expressed in terms of /j. 

The target nucleon is characterised by certain quantum numbers such as 
isospin and baryon number. These quantum numbers are carried by the con- 
stituents. For example, the net number of up quarks in a proton is two, hence 
Jq 1 dz(f u (z) — fu(z)) = 2. By forming appropriate combinations of scattering 
cross-sections, quantities can be measured that correspond to these conserved 
quantum numbers and hence have simple values in the Naive Parton Model. 
These quantities are referred to as Parton Model sum rules. They can then be 
compared with experiment and the fundamental properties of theory tested. 
This Naive Parton Model is subject to corrections in the full theory of strong 
interactions (Quantum Chromodynamics or QCD |3], |J, for recent reviews 
see H [| [7| ). These corrections fall into two types; those that are strongly 
suppressed at high energy (higher twist corrections) and those that vanish 
only logarithmically with the momentum transfer. The latter are fully cal- 
culable in terms of the coupling constant a s of QCD. Comparison of the sum 
rules with these QCD expectations then provides a powerful test of QCD 
and enables a s to be measured. 

In the remainder of this article we will discuss these sum rules. We show 
their values in the Naive Parton Model and the corrections from QCD and 
finally compare these predicted values with experiment. 
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Figure 1: Diagram illustrating the deep inelastic scattering process electron 
+ proton — ► electron + anything. 

1.1 Deep Inelastic Scattering 

The kinematics of lepton-nucleon scattering are shown in Figure 1. The 
scattering of an unpolarized charged lepton (an electron or muon) or neutrino 
of momentum k off an unpolarized nucleon of mass M, momentum P results 
in a final state with a lepton of momentum k! and a nuclear fragment. The 
inclusive cross-section for this process can be described by two Lorentz scalars 
which can be taken to be Q 2 — — {k' — k) 2 and v = P ■ (k — k'). Convenient 
dimensionless variables are x = Q 2 /2v and y = v j P ■ k both of which range 
from zero to one. The latter is the fractional energy loss of the lepton in the 
rest frame of the nucleon. Two types of scattering are important, neutral and 
charged current scattering. The former case describes the process fiN — > \iX 
or eiV — > eX. For neutral current scattering and Q 2 « M§ (all of the data 
discussed in this article satisfy this requirement), the cross-section can be 
expressed as 

dxdy sx 2 y 2 ^ 
(xy 2 F? N (x, Q 2 ) +(l-y- M 2 x 2 y 2 /Q 2 )Ff{x, Q 2 )) 

where s = (P + k) 2 and F 1 , F 2 are arbitrary functions called structure func- 
tions. This is the most general form that this cross-section can take con- 
sistent with Lorentz invariance and parity conservation provided that terms 
that are proportional to quark masses (actually m 2 /Q 2 ) are neglected. This 
is an excellent approximation since the nucleon consists mainly of up, down 
and strange quarks that have very small masses (m u ~ m d ~ few MeV and 
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m s ~ 100 MeV ||). For charged current scattering, the presence of parity 
violation which proceeds via the exchange of a W boson between the lepton 
and the target nucleon, requires the appearance of another function F 3 . 

dxdy" 2u{Q 2 + M 2 ,) 2 ^ 1 [X > ^ } 

+(i-y-x 2 y 2 M 2 /Q 2 )F» N (x,Q 2 ) (2) 

For V N scattering the sign of the last (xF 3 ) term is reversed. The above 
formulae assume that the target nucleon and charged lepton are unpolarized 
(the neutrino is always polarized). In the case of eN scattering with polarized 
electron and nucleon a parity conserving asymmetry can be formed. 

do- eN dcj eN 

a U y) = a(Ja P ( 3 ) 

dxdy dxdy 

where the subscript p (ap) refers to the state where the nucleon spin is parallel 
(anti-parallel) to its direction of motion in the center of mass frame of the 
lepton-nucleon system. In both cases the lepton has its spin aligned along 
its direction of motion. Then 

a(x, y) = ^#((1 - 2/y 2 + 2x 2 y 2 M 2 /g 2 )G 1 (x, Q 2 ) 

+Ax 2 M 2 G 2 (x,Q 2 )/Q 2 ) 

defines the spin dependent structure functions G± and G 2 - 

In the rest frame of the target nucleon, the various kinematic quantities 
are related to the energy E (E') of the incoming (outgoing) lepton and the 
scattering angle 9 of the lepton by y = u/ME = (E — E') / E and sin 2 (#/2) = 
Q 2 j (4EE'). By measuring the incoming and outgoing lepton energy and the 
scattering angle the structure functions can be determined. 

1.2 Quark Parton Model 

In the Quark Parton Model, lepton-nucleon scattering is described by the 
scattering of a lepton off the partonic (i.e. quark and gluon) constituents of 
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the nucleoli. The nucleon structure is described in terms of the parton distri- 
bution functions fi(z). We shall often use the following symbols to simplify 
the notation :- fi(z) = i(z) etc. so that the up quark distribution in a proton 
is u(z) and the anti-down quark distribution is d(z). These distributions will 
always refer to a proton. When a neutron target is involved we use isospin 
symmetry to relate the neutron distributions to those of the proton so that 
the up quark distribution in a neutron is the down quark distribution in the 
proton and vice-versa, while the gluon (g(z)) and the strange, charm, bot- 
tom and top distributions (s(z),c(z),b(z),t(z)) are the same in proton and 
neutron. The distribution of these heavier quarks is smaller and in the case 
of top and to a lesser extent bottom totally negligible. The charm quark 
is troublesome since its mass cannot be neglected in experiments that have 
Q 2 ~ M 2 ~ (1.5 GeV) 2 . This effect is most important in neutrino scattering 
where the process u + s — >/i~ + c is a significant part of the cross-section. 

Since the quarks carry the quantum numbers of the nucleon these distri- 
bution functions satisfy certain constraints. For example, the electric charge 
of the proton yields 

1 = jf 1 dz{\{u{z) - u(z)) - l -{d{z) - d{z)) - l -{s{z) - *(*))) (5) 

and the zero net strangeness of the proton gives 

= f 1 dz(s{z) -s(z)) (6) 

J 

Note that this does not imply that s(z) = s(z). Momentum conservation 
in the lepton-parton scattering process implies that the parton momentum 
fraction z is identified with the kinematic variable x in the Naive Parton 
Model. The structure functions F\, F 2 and F3 are given in terms of the 
parton distribution functions. In particular 

F 2 (x,Q 2 yv = xY,q-Mx) (7) 

i 

where g« is the charge of the type % parton. The Naive Parton Model therefore 
predicts that there is no Q 2 dependence in Q 2 ). The partons that couple 
to the photon and W boson (quarks) have spin 1/2 and hence F 2 (x, Q 2 ) = 
2xFi(x,Q 2 ), the Callan Gross relationship ||. 



4 



The relationships between the structure functions for neutrino scattering 



and the parton distributions are complicated by the Kobayashi-Maskawa flit) 
mixing matrix, tyy which determines the relative strength of the coupling of 
a quark pair to the W boson :- W^qff^l — JsjqjVij. In order to simplify 
the results that follow, the equations are written in the approximation that 
the mixing matrix is a diagonal unit matrix. The mixing can be added 
by replacing the quark distributions by appropriate linear combinations, for 
example d(z) — > (Vd u + Vdc + Vdt)d(z). Neutrino scattering proceeds off up 
and anti- down type quarks viz vu — > fi + d and vd — > (i + u leading to the 
following relations (neglecting top and bottom quark contributions) 

F^(x, Q 2 ) = 2{d{x) - u{x) - c{x) + s{x)) (8) 

and 

Fz p {x, Q 2 ) = 2x(u{x) + c(x) + d{x) + s(x)) (9) 

The other relations can be trivially obtained from these. The charm quark 
contributions will not be written explicitly in the following. 

The spin structure of the nucleon is probed in polarized scattering. One 
can define Afi(z) as the difference between the parton distributions for par- 
ton of type i with spins (helicity) parallel and anti parallel to the nucleon's 
spin. The unpolarized distributions introduced above are then the sum of 
these two states. Define the quantities clq, a 3 and as by 



a 

«3 



J dx(Au + Au + Ad + Ad + As + As 
J dx(Au + Au - Ad - Ad 



'0 

1 r 1 

dx[Au + Au + Ad + Ad-2As-2As 



(10) 



These quantities are related to the matrix elements of the axial vector current 
between nucleon states. The matrix elements involving changes of flavor can 
be determined from weak decays. (13 is determined from neutron /3-decay: 
a 3 — 9a = -q^ = 1.2573±0.0028. a§ and 03 can be constrained from the weak 
decay constants of (E, A and S) hyperons 



n . Assuming SU(3)p symmetry 
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for the octet axial vector currents this gives F = (V3a 8 +a 3 )/4 = 0.459±0.008 
and D = y/3(y/3a 3 - a 8 )/4 = 0.798 ± 0.008 implying a 8 = 0.33 ± 0.02 g|. 
Data from up and Vp elastic scattering |T3f provide information on the matrix 



element of the flavor singlet current ao to be determined [14], These data 
do not directly measure the static quantities and a form factor behaviour 
must be assumed ||16|| . In addition the experimental errors are large. The 



data can be interpreted as Jq dx( As + As) = AS = —0.15 ± 0.09 [[L4|] which 



(using the above value of a§) implies ao = 0.12 ± 0.27. If all of the nucleon's 
spin is carried by quark spin and not by gluons or orbital angular momentum 
one expects that = 1. In the Parton Model the spin structure functions 
are 

G 2 (x,Q 2 ) =0 

(11) 

G?(x,Q*)= iEi^A/ 4 (x) 

1.3 Sum Rules 

In this section we list the various sum rules and their values in the Naive 
Parton Model. These rules are all derived from inclusive quantities that 
have a simple interpretation in this model. We will refer to these sum rules 
by a name that relates to this simple interpretation. We will also indicate the 
more familiar names by which they are sometimes referred in the literature. 



The baryon (Gross Llewellyn- Smith) sum rule |T7| uses the av- 
erage of F3 measured on a proton and a neutron. Since neutrino experiments 
are often performed on heavy nuclear targets (such as iron) which have an al- 
most equal number of protons and neutrons, the quantity is readily measured. 
The sum rule measures the sum of the baryon number (B) and strangeness 
(5*) of the nucleon. 



S GLS = f 1 dxFf = dx\ (F£ p + F> 
Jo Jo 2 

dx(u(x) + d(x) — u(x) — d(x) + 2(s(x) — c{x))^j (12) 



The Isospin (Adler) sum rule Jl8[ measures the isospin of the target and 



6 



depends on the difference in F2 measured in neutrino scattering off proton 
and neutron targets. 



S A = [ — (FP - F» p ) 
Jo x 

= 2 J dx(u(x)-d(x)-u(x) + d(x)) ( 13 ) 
= 4/3 = 2 

A similar sum rule can be formed in electron scattering. We can introduce 
"valence" distributions defined by u v (z) = u(z)—u(z) and d v (z) = d(z)—d(z). 
Note that, since the net number of up (down) quarks in a proton is 2 (1), 
these satisfy / dzu v (z) = 2 and / dzd v (z) = 1. Then 

S G = f 1 — (F? - fd 
Jo x 

If 1 / _ - \ 

dx[u{x) — d(x) + u(x) — d(x)j 



3 Jo 

I rl , v 2 f 1 (_ \ 

- / dx[u v {x) — d v (x) ) H — / dxiuix) — d(x) ) 
3 Jo v ' 3 jo ' 

1 



(14) 



The last step follows if / dzd(z) = J dzu(z). However, there is no fundamen- 
tal reason for this simple assumption to be valid and hence the Valence 
Isospin (Gottfried) sum rule |L9| is on much weaker ground than the 
previous one. 

The total momentum carried by all of the proton's constituents is con- 
strained to add up to that of the proton. Hence the Momentum sum rule. 



S'mom = Jq 1 dx x (u(x) + d(x) + s(x) + u(x) + d(x) + s(x)j 
= 1 — Jq 1 dx xg(x) 



(15) 



This sum rule cannot be tested directly since the gluon distribution function 
does not appear in the structure functions. Rather, the left hand side of this 
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equation can be measured and the rule used to infer something about the 
gluon distribution. 

The Second Isospin (Unpolarized Bjorken) sum rule [EU1 is 



^unpol. = £ dx ^ p un _ pvp^ 

dx[u(x) — d(x) — u(x) + d(x)J (16) 



It is equivalent to the Isospin SUM rule in the Naive Parton Model where 
F 2 (x,Q 2 ) = 2xFi(x,Q 2 ). This equivalence is broken once QCD corrections 
are computed. 

The remaining sum rules depend on the spin structure function G\(x, Q 2 ) 
and therefore on the polarized distribution functions. The integrals are re- 
lated to the quantities ao, a 3 and ag by 

The Polarized Isospin (Bjorken Spin) sum rule pUj : 

S Bi = C dx{G{-G n 1 ) = \a z (18) 
Jo o 

If we assume that the strange quarks are unpolarized As = — > ao = V3a$ 
then we have the Spin (Ellis Jaffe) sum rules p2| |3 



(19) 



The Parton Model prediction G 2 (x,Q ) = leads to a trivial prediction for 



the G 2 (Burkardt-Cottingham) sum rule p4]. 



S BC = [ dxG 2 (x) = (20) 
Jo 
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2 Sum Rules in QCD 



In the full theory of strong interactions (QCD), the Naive Parton Model 
and its expectations for the values of the sum rules are modified. These 
modifications are of two types. At high energy (large momentum transfers), 
the coupling strength of QCD becomes small and perturbation theory can 
be used @, [§]. In this regime, corrections to the sum rules can be expressed 
as a power series expansion in the strong coupling constant a s (Q 2 ). At lower 
values of Q 2 non-perturbative corrections enter which can be expressed as 
power series in 1/Q 2 . Unlike the perturbative corrections, these cannot be 
calculated at present. In some cases, corrections to different processes can 
be related to each other and experimental results may be used to determine 
the effect of these corrections on the sum rules. This section analyses both 
the perturbative and non-perturbative corrections. 

In Section 2.1 perturbative corrections are given first in the framework 
of the QCD improved Parton Model. This approach allows to give a very 
appealing and intuitive description of the basic ideas of factorization and 
Q 2 evolution of structure functions. We restrict our discussion in this para- 
graph to leading order corrections and will see that in the leading logarithmic 
approximation all sum rules remain valid. 

For the discussion of higher order perturbative corrections it is conve- 
nient to employ the framework of the operator product expansion, since the 
structure of the corrections becomes most transparent in this more formal 
approach. It leads of course to the same results as one would get with the 
QCD Parton Model. In fact, the connections between both descriptions will 
be pointed out wherever possible. 

The operator product expansion has the further advantage that non-per- 
turbative effects can easily be incorporated. Power corrections of higher twist 
are studied in Section 2.2. 

2.1 Perturbative QCD Corrections 
QCD Parton Model 

In the previous chapter the relation between the structure functions and 
quark distributions was given and it was stated that the Naive Parton Model 
predicts the Q 2 -independence of the structure functions F(x, Q 2 ). Violations 
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of this scaling behaviour were observed experimentally and may be explained 
theoretically due to strong interactions. The QCD improved Parton Model 
gives a simple and quantitative description of these effects and introduces 
the correct Q 2 dependence into the parton distribution functions. 

The QCD generalization of Eqn. (^) is provided by the factorization 
theorem for deep inelastic scattering (see |25|, |26j 

^ NS (z,Q 2 ) = f^cJ^^a^))^,^,^). (21) 
Jo y \y n 2 fx 2 ) 

Here serves as a generic notation for 2Fi,F 2 /x and F 3 . For polarized 
structure functions T\. = 2Gi, G 2 one should replace fk by Af k . For simplic- 
ity we consider the nonsinglet combination of structure functions like ep — en. 
Accordingly f^ s are the appropriate combinations of parton densities, e.g. 
for = F 2 NS /x one has / 2 NS = J2iQiff S - In Eqn.(gTj) the factorization of 
high momentum (short distance) and low momentum (long distance) effects 
is expressed. The former are described by the coefficient function and cal- 
culable in perturbation theory. As a characteristic feature for perturbative 
computations one finds a dependence on the renormalization scale /i 2 . Long 
distance contributions cannot be calculated by present theoretical methods 
available in QCD and are absorbed in the parton distribution functions. The 
separation between the low and high momentum regime calls for another 
scale, the factorization scale fif. In the following we shall choose fij = /i 2 . 

The coefficient functions of Eqn.([H]) were calculated in the leading log- 
arithmic approximation in [ 2"7f , to order a s in |p8 , 



and order a? in 



[ PH| , |3T] , |32] , 33, 33]. In leading order QCD they have the following form 
(in this and all subsequent equations a s without an argument is understood 
to mean a s {Q 2 )) 



G 



NS 
k 




X 



Q 2 



P qq \-)\n^ + R 



->NS 



(22) 



The well known splitting function P qq (z) = C F [(l + z 2 )/(l-z) + + (3/2)5(1 - 
z)] measures the variation with Q 2 of the probability of finding a quark 
inside a quark with a fraction z — x/y of its momentum y ]27[|. This leading 
logarithmic term is universal for all structure functions and can be absorbed 
into newly defined, Q 2 -dependent parton distribution functions 



/r(*,Q 2 )=/ fe NS (*)+ 



2^ 



f A NS (0 



P qq \n% + R™ 
/i 



fc,abs 



(23) 
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The particular way how R^ s is split into an absorbed part -R^lbs an d a 
remaining part R^ em is a matter of convention and specifies the so-called 
factorization scheme. Two popular choices are the DIS scheme and the MS 
factorization scheme (for details see e.g. [p6fl). With the Q 2 dependent quark 
distributions the structure functions can be rewritten in the following form 
if terms of order a 2 are neglected: 



^ NS (*,Q 2 ) = —15(1--] + y-Rk% m ) fk S (y> Q 2 ) 

Jo y V V V) 27r / (24) 
L -^fr(x,Q 2 )+0(a s ) 

It can be seen in the last step that in the leading logarithmic approxima- 
tion (LLA) the relations between structure functions and parton distribution 
functions remain unchanged except that the parton densities now depend on 
Q 2 . With the same modification all Parton Model sum rules of the previous 
chapter remain valid in this approximation. 

The Q 2 dependence of the distribution functions and hence the struc- 
ture functions and sum rules is most readily expressed by the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) |35| , |27|1 evolution equation which 
to leading order in QCD has the following form 

Q2dJl ^) ^fa^r^+toWrf)*. (25) 

If one is not restricted to flavour non-singlet combinations the other function 
Pqgiy) comes into play due to the probability of finding a quark inside a 
gluon. Both splitting functions are determined by perturbative QCD and can 
be written as an expansion in a s . This equation can be used to determine the 
perturbative QCD corrections to the various sum rules. It should be noted 
that the DGLAP equation contains more information about the behaviour of 
the structure functions than does the set of sum rules. However, the higher 
order QCD corrections to the sum rules are easier to compute than those to 
the DGLAP equation. Hence while the DGLAP evolution equation is only 
known to order a 2 0, ^ [H| §§, @, g§ g§ g5|, the corrections 



that we discuss below to some of the sum rules are known to order a^. 
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Operator Product Expansion Approach 



Modifications of the sum rules due to higher order QCD corrections were 
indicated in the previous section for the QCD improved Parton Model. In 
this section we explicitly discuss those corrections. Their structure becomes 
particularly transparent in the framework of the operator product expansion 
[ fl6]| . This approach will also prove to be useful in the following section for the 
discussion of non-perturbative effects. The Mellin moments of the structure 
functions are expanded in a form [|7|] where the short distance and long 



distance contributions are factorized in a similar fashion as in Eqn.( 



o 




dx x^F k (x, Q 2 ) = V Cif n i=L aB tf) Agtf). ( 26 ) 



The expansion is expressed in terms of reduced matrix elements A^f of opera- 
tors renormalized at scale \i which describe the long distance effects and have 
to be determined from experiments and coefficient functions C^ T n that de- 
scribe the short distance effects and can be calculated perturbatively. The op- 
erators are characterized by their quantum numbers which can be broken into 
two types. The label r refers to the twist of the operator and i refers to the 
flavor quantum numbers such as isospin. The "twist" of an operator is defined 
by its dimension do minus its spin n. Since G % kn (Q 2 /fx 2 , a s ) ~ (/i 2 / 'Q 2 )( T ~ 2 ^ 2 , 
operators of lowest twist (r = 2) dominate in the large Q 2 limit. We shall 
postpone the discussion of higher twist operators to the next section and 
omit the twist label when considering the leading, twist-2, terms. 

If one considers for simplicity the moments of flavour non- singlet struc- 
ture functions, one can readily see the similarity between the approaches of 
the QCD Parton Model and the operator product expansion. Taking the mo- 
ments on both sides of Eqn.(^Tj) shows that the Parton Model analogues of 
the coefficient functions and operator matrix elements of Eqn. ( p6|) are given 
by the following moments 



Cjg^a.^ 2 )) <=► f Q dxx^Cr[x^,a s {^ 



(27) 

Jo 

The Q 2 behaviour of the coefficient functions is governed by their renor- 



malization group equation IS , which follows from the fact that the LHS of 
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Eqn. (p6|) as a measurable quantity is independent of /i 2 . The non-singlet op- 
erators are renormalized multiplicatively with renormalization constant ■ 
The renormalization group equation then has the following simple form 



\~~2J A(fJ, 2 ) 



o. 



(28) 



where A = a s /ir. The beta-function and the anomalous dimension are de- 
fined by 

d A / o\ Ng o d 



d/r 



log 3? 



NS 



and are given by the expansions 



(29) 



{3(A) = -p A 2 - P,A 3 - p 2 A A + ... 
7f (A) = 7f (0 M + 7 f (1) A 2 + 7f (2) A 3 + • • • • 



(30) 



The solution of Eqn.(|28|) is 



C% iX> A{f)\ = C™ (l, A(Q 2 )) exp I f_ 



(31) 



where the formula for the effective coupling constant d s and the coefficients 
of the anomalous dimension and the beta-function are listed in the appendix. 

From Eqn. (|3l|) it becomes apparent that QCD corrections in are 
twofold. They arise from the expansion of 



cz (i, Am) = i + BZ (1) Am + b^a\q* 



(32) 



as well as from the expansion of the exponential on the RHS of Eqn. (|31|) . 
The Q 2 dependent part of that term is 



*m t ns (2;) 

exp ( / dx— 

■U + A{Q 2 ) 



A(Q* 



NS(0) 
[n 



7„ NS(1) A7„ NS(0) 



'o 



(33) 



+ A 2 (Q%..} 
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The leading term 7^ S( ^ of the anomalous dimension is independent of the 
renormalization scheme. This is no longer the case for the higher order 
terms ^ s ^\i > 0. However, in expressions for physical quantities they are 
associated with the scheme dependent coefficients B^' n in such a way that 
the final answer is renormalization scheme invariant. 

The /i 2 dependent part of the exponential in Eqn. ([IT]) may be combined 
with the reduced matrix element A^ s (fi 2 ) to form the renormalization scale 
invariant expression 

A Ns,inv = exp f f A ^ 2) d x^M\ Af (/i 2 ). (34) 



(3{x) 

In the case of the sum rules all relevant operators have 7„=i = 0. This can 
be observed from the relation Jq dx x n ~ x P qq (x) = 7^ s (°)/2, which vanishes 
for n = 1 as as consequence of fermion number conservation. 

In view of Eqns. (|32| , |33|) the RHS of Eqn. fl2"6|) approaches a constant value 
as Q 2 — > oo that is basically given by t4^ s ' 1i1v and may be identified with the 
corresponding expression obtained in the Naive Parton Model. As will be 
seen below the situation becomes even simpler, when the operators under 
consideration are conserved currents. In this case the anomalous dimensions 
vanish and the QCD corrections are already completely determined through 
Eqn.([3^). Since these operators are not renormalized, their reduced matrix 
elements A^ s are independent of \j? . 

Analogous relations to Eqns. ( p8| , |31| , |34]) hold for singlet combinations of 
structure functions. They are more complex than for the nonsinglet case, 
since mixing of different operators with the same quantum numbers may 
occur under renormalization, leading to an anomalous dimension matrix. 

Let us illustrate the above discussion in an example and consider the 
moments of the structure functions Fi and G\ of deep inelastic electron- 
nucleon scattering 

f 1 dx x- 2 Fr en = E CknV? S n = 2, 4, 6 . . . 

ri 1 (35) 

I dx x"- 1 Gf =-J2 Ch un af n = 1, 3, 5 . . . 



which leads to the Valence Isospin sum rule and the Spin sum rules 
respectively. Depending on the crossing properties of the structure functions 



14 



under fi «-> v y x «-> — x, only operators with definite spin signatures are 
relevant in the expansions of Eqn.(|35]). Spin-even operators contribute to 
the moments of the combination p^ p ~ en anc l spin-odd operators to those 
of G\ N . For later use we note that F% l £~ vp , p" n+y P a i so represent spin-odd 
combinations. 

It is obvious that the Spin sum rules are immediately obtained from 
the first moment n — 1 of G\ in Eqn. (|35|) . In this case the corresponding 
operators are the flavour nonsinglet and singlet axial vector currents 

= <H*7s ^yj V> ; J 5° M = ^75^ (36) 

of the SU (3)f symmetry group (j = 1, . . . , 8) (the Xj are the Gell-Mann ma- 
trices). The af are given by the matrix elements of these operators between 
the states of the nucleon N = p,n with momentum P^, spin and mass M 

(P,S\J^\P,S) = MafS fl ; ; (P, S\ J°JP, S) = 2Ma$S„ (37) 

The nonsinglet axial vector currents are conserved in the massless quark 
limit. According to the discussion above the are independent of fi 2 and 
the corresponding anomalous dimension vanishes A7„j! 1 = 0. (We use the 
notation A7„ for polarized scattering in distinction to 7 n for unpolarized 
scattering.) The scale dependence of a (fi 2 ) on the other hand reflects the 
fact that the singlet axial vector current is not conserved due to the ax- 
ial anomaly |49|, [5(| and therefore has a nonvanishing anomalous dimension 
A7^ =1 7^ 0. The polarized anomalous dimension A^ s ^ s for arbitrary n 



was calculated in leading order in [|5T], 0, [27J and next-to-leading order in 
3^, [34], j54|. The next-next-to-leading order result &ln,qq can be found in 



3^, |56[. The nonleading results were calculated using the MS scheme, the 
standard modification of the Minimal Subtraction scheme (A notewor- 
thy feature of this particular renormalization scheme is the separate gauge 
invariance of both the anomalous dimensions and the coefficient functions 
[01 .) Some of them are listed in the appendix. The Naive Parton Model 
expressions for the quantities ag nv , a 3 and a 8 as well as their values were 
already presented in Eqn. (|To|). 

The situation for unpolarized lepton-nucleon scattering is not as straight- 
forward as for the polarized case. Since only spin-even operators contribute 
to the operator product expansion for the moments of F 2 in Eqn.fl3"5|), the 
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Valence Isospin sum rule cannot be obtained simply as the special 
case n = 1 of the first moment. Nevertheless, the information that is con- 
tained in Eqn.(^) can be used to derive QCD corrections to the Valence 
Isospin sum rule. The quantities v^ s are defined as the reduced matrix 

elements of operators of the form of 8 " 1 '-"" = {^7 W D>» ..D^(X j /2)ip} 

with n = 2, 4, 6, . . . and where S indicates symmetrization of the indices. The 
perturbative expansions for the corresponding anomalous dimension 7^ and 
the coefficient function C|J n (l, A(Q 2 ) are a priori meaningful only for even n. 
However, the QCD Parton Model not only reproduces the results of the op- 
erator product expansion, but also provides an answer for moments with odd 
n, for which no operators are available. QCD corrections to the Valence 
Isospin sum rule may therefore be obtained by analytically continuing the 
results valid for even n to the formally forbidden values of odd n. Similarly 
a continuation for spin odd combinations of structure functions to values of 
even n can be made. Consequently the generalization of the anomalous di- 
mensions 7n=evcn/odd are denoted as 7*, now valid for all n. In leading order 
the nonsinglet and singlet anomalous dimensions 

7 NS/S(0) we 

obtained in 

IPI |60| , p7[| . The nonsinglet (singlet) next-to-leading order MS result 7^ IS( - 1 ^ ± 
(7^) was calculated in pi (|B2L SSL B3) and simplified in §3 (p§). Fur- 



ln,%3 

thermore, the anomalous dimensions are given by the moments of the split 
ting functions P^p , which were directly computed in |37|, f4T| 



Three loop nonsinglet anomalous dimensions 7^ 2 ' ) are given in for 
even n > 2. We note in passing that spin-even and spin-odd combinations 
of structure functions correspond to combinations qj ±qf and Aqf =p Aqf of 
parton densities (see e.g. 0, IBB]). 

We have demonstrated that QCD corrections to the Valence Isospin 
SUM rule can be derived from the operator product expansion for the mo- 
ments of F2, even though no operators exist for the first moment. That there 
is no corresponding operator for n = 1 is reflected in the fact that no reliable 
numerical value for Sq can be given. Indeed, the Parton Model prediction 
for the Valence Isospin sum rule was based on an additional assumption 
about the quark sea which has no solid theoretical justification. 

Within the same flavour octet the Q 2 -evolution for all nonsinglet op- 
erators in Eqn.(^Bj) is given by a common coefficient function defined by 
Cl,n = ^fcC^n; i = 1..8. Similarly, for the first moment one may denote the 
singlet coefficient function by C^ n=l = <5°Cf. The constants 8\ are combina- 
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tions of quark charges or weak couplings. In our example of the neutral cur- 
rent (J e ™ = EfQ fifing) Process one has 5 3 k = 1/6, 5 8 k = 1/(6^3),^ = 2/9 

G\ . Next-to-leading order coefficient 



and all other 51 



for k 



T?eN 



functions for the structure functions Fi, F 2 , F 3 of electron and neutrino scat- 
tering off nucleons are given for the MS scheme in [|9|, |62], [70], [Ti]] and for 
other schemes in J72|, [73|, [F4[] , |75], [76], |28|, [f7], [ft| . For the various sum rules 



only the numerical values of the corrections will be given in this section. 
The corresponding analytic expressions for the coefficient functions and the 
anomalous dimensions can be found in the appendix. 
For the Valence Isospin sum rule 



S G 4c & --,^,Q 2 ,,ex P (/ J, ° 2, ^2|W- 



(38) 



it was found |]69| that B^n=i = vanishes. The QCD corrections are there- 



NS(1)- 



fore due to the nonvanishing anomalous dimension 7 n= i 



S, 



G 



1 + 



1 13 + 8C(3) - 2vr 2 a 



33 -2n 



f 



7T 



0.036\ f<\ 
0.0387 



7T 



(39) 



Here the upper and lower coefficients refer to rif = 3 and 4 respectively. 
For polarized electron-nucleon scattering the moments read 



dxG 1 



ep(en) 



; I) ^ = q s 1 (i 1 ^))i(±a 3 + 1 



+67^(1, A(Q 2 ))exp 



-Am A 7 s =1 



I AW) 



Cg 1 S i(1^W 2 ))^(±«3 + ^ 



A(Q 2 ))exp 




"A(Q 2 ) A 7 



(40) 



In the nonsinglet part of this equation the exponential does not contribute 
because of the vanishing anomalous dimension A7^ 1 = 0. In the singlet 
part the renormalization scale dependent piece of the exponential may be 
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combined with the operator matrix element to form the scale independent 
quantity a™. The corrections to the nonsinglet coefficient function Cq S u=1 

2 §g, §TJ, and al @ and for the singlet 
D, B3I. The QCD cor- 



were calculated to orders a* |T7D 



coefficient function Cq n=1 to order a s [^3 



and 



rections for the Polarized Isospin SUM rule are completely determined 
by the nonsinglet coefficient function CW n=1 . The numerical result reads 



S 



Bj 



6 



1 _ ^1 _ 
7T 



3.583\ 
3.250/ 



7T 



/20.215 
U3.850 



7T 



(41) 



The Spin sum rules are obtained by inserting the coefficients and 
anomalous dimension of the appendix and read 



a 



p(n) 
EJ 



1 

12 



3.583 
3.250 



[ ±a3 + 7S a \ 

~ 12 [ 7T 

±a 3 



1 



+ gao(Q 2 ) 



1 



7T 



/20.215 
U3.850 
1.096 \ 
-0.067/ 



3.583 
3.250 



as + 



9 



1 - 



a s \ 2 _ /20.215 
IFJ ~ U3.850 
0.333\ a s _ 
0.040/ ~ ~ 



7T 

71 

0.550 \ /a t 
-1.082/ \ tt 



(42) 



where in the first step the choice fi 2 = Q 2 was made. 

The validity of the G2 SUM rule can be derived from the nonexistence of 

Mr 



leading twist operators for n 



8q| in the operator product expansion 
of the moments of G^- A continuation from higher n to n = 1, similar to the 
case of the Valence Isospin sum rule, also leads to a vanishing result due 
to a kinematical factor (n — l)/2n in the expansion (see [0). This result 



is confirmed by explicit calculations in |86], B7|| , despite an earlier contrary 

(43) 



claim B8 



S 



BC 







We now turn to the situation for neutrino or antineutrino scattering on 
protons. The combination F^ p + F% n in the BARYON SUM RULE transforms 
as a flavour singlet. In view of our earlier discussion one could expect a 
more complex renormalization group equation than Eqn. (p8|) for CpJ^ =1 . 
However, in this case the gluon field operator transforms differently under 
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charge conjugation than the quark field singlet operator and therefore no 
mixing occurs. Since the anomalous dimension vanishes, the QCD corrections 



are due to the coefficient function Cp s J |613|, |8~2|, |3~3~ 



MO 



S, 



GLS 



7T 



3.583 
3.250 



2 / 18.976 
~~ U2.198 



(44) 



It was pointed out in [ISO, ISTl, [3~1] that the corrections up to second order 



a: 



for the Baryon sum rule are identically the same as for the Polar- 
ized Isospin SUM rule. As may be seen from the explicit formulae in 
the appendix, corrections for both sum rules differ at order due to terms 
~ d abc d a bc, representing diagrams with an internal fermion loop and a purely 
gluonic intermediate state. 

For the Isospin SUM rule both the anomalous dimension and the cor- 



rections to the coefficient function B 



(vn—vp) 
F 2 ,n=l 



up to order a 2 s were found to 



vanish [BE, 32]. Therefore the Isospin sum rule is not violated: 



Sa = 2 



(45) 



Since this sum rule was derived by the use of current algebra methods [|ilf|, 
it is expected to hold true in QCD for all orders in a s . 

The Second Isospin sum rule is also characterized by its vanishing 
anomalous dimension t^Ej = 0, but has a nonzero coefficient function [69 
|89|, |9(| . As a result one has 



riunpol 



1 - 0.667- 



7T 



2.944 
2.648 



/18.596 
U3.381 



(46) 



Combining the Isospin sum rule and the Second Isospin sum rule 
leads to the corrections for the first moment of the Callan-Gross Rela- 
tion for neutrino scattering. 



1 n T 

""^ T-ivn—vp 

b L 

X 



f-(F.. 
Jo X 



vn—vp 
2 



2xF{ 



vn—vp\ 



2 cjiunpol. 



(47) 



The Callan Gross Relation for electron-nucleon scattering is violated in 
QCD as well. The coefficient functions of the longitudinal structure function 
were calculated to leading order in |72|, [75], |9l], [73], |74]|. The results for next-to- 
leading order corrections turned out to be very controversial. Nonsinglet or 
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singlet coefficient functions Cl{x, Q 2 ) and their moments CL, n {Q 2 ) were given 

) and PI |97|, I 



51 30, 31 



(see also 



in order a 2 in ]92], |93| , 
respectively. There is complete agreement between 



99, 100 



and pq 



Third order 0(a^) corrections to the second and higher (even) moments of 
the nonsinglet longitudinal structure function can be found in p6|, ^7|. The 



corrections to the Callan Gross Relation enter the cross-section ratio 
R = <jl/(Tt which is measurable in experiments. 

There have been attempts to estimate even higher order corrections. Ref- 
|101|1 estimates the 0{a i s ) coefficients for the Second Isospin sum 



erence 



RULE and the nonsinglet part of the Spin sum rules. Neglecting light-by- 
light diagrams, this estimate is also used for the Baryon sum rule. These 
results are in reasonable agreement with another estimate by [ |102 ] using Pade 
approximation. Reference ||103|| gives the 0{a z s ) estimate for the singlet part 
of the Spin sum rules. 



2.2 Higher Twist Effects 

So far only leading twist corrections to the sum rules have been discussed and 
nucleon mass effects were neglected. Contributions of higher twist operators 
in the operator product expansion of the structure function moments are 
suppressed by powers of l/Q 2 - These so-called dynamical power corrections 
are of non-perturbative nature and may become important at low Q 2 . They 
are difficult to estimate, because at present no method is available for reliably 
calculating in the non-perturbative regime. Estimates of higher twist effects 
are therefore usually accompanied by large uncertainties. 

Another class of power corrections are of purely kinematical origin and 
arise in case that the nucleon mass is not neglected. Such target mass effects 
are suppressed by powers of M 2 /Q 2 and may also become relevant for low 
values of Q 2 . They can be taken into account to all orders in (M 2 /Q 2 ) n 
by replacing the structure function moments through so-called Nachtmann 
moments | 104|| . 



A standard approach ||105| , |10(j|| to calculate higher twist contributions 



is based on QCD sum rules ||107|| . This method parameterizes nonpertur- 
bative effects in terms of quark and gluon condensates. These condensates 
are matched with a hadronic representation obtained phenomenologically 
through a dispersion integral with a spectral density given by resonance and 
continuum contributions that can be extracted from data. The matrix ele- 
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ments of the higher twist operators are therefore expressed through conden- 



sates, the values of which, however, are only known to 10-20% [107 



Power corrections ~ l/Q 2 to the first moments of the polarized structure 
function Gi are given by [|108|| (see also [jS5fl ) 



where 



((O)) 



9 



M 2 

((U))-—((V)) 



and {{U)),((V)) are defined by 



(48) 



(49) 



(P,S\U,\P,S) = 2MS,((U)) 
(P, S\V^\P, S) = 2M((V)) {(S„P V - S„PjP a )} S{u>a} 
as the reduced elements of the operators 



(50) 



f/S/NS 



V UT = 9s {uG^u ± dG^d} a 



(51) 



for the proton (u <-> d leads to the corresponding neutron operators). Here 
the strange quark contribution is neglected and 



G^ u — -e 



(52) 



After some errors of ||106|| were discovered in [109] , the analysis of higher 
twist effects in was updated (see [§§]). 



0D EJ 



0D EJ 



(0.02±0.013)GeV 2 
Q 2 

(0.005 ± 0.003)GeV 2 
Q 2 



(53) 
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These values agree approximately with another analysis employing the QCD 
sum rule approach ||1 10|, llT|. They are, however, considerably smaller than 



those obtained in ||112| , |113j , |114| | which rely on the relation between the 
first moment of the polarized structure function G^ n and the Gerasimov— 
Drell-Hearn sum rule ||115| , |116|| for forward scattering of real photon off 
nucleons. A re-examination of the QCD sum rule method in ||117|| arrives at 



-(0.029 GeV 2 )/Q 2 and SS% 3 = -(0.002 GeV 2 )/Q 2 with 



the result SS^j 

uncertainties of about 20%. A recent paper [ |118|| calculated the one-loop 
anomalous dimension for the r = 4 operators that allow the logarithmic 
corrections to the 1/Q 2 term to be computed. 

Power corrections for the Polarized Isospin SUM rule follow imme- 
diately from Eqn.( 



6S, 



Bj 



The first term was estimated in [fl0§ to be -(0.09 ± 0.06)GeV 2 /(6Q 2 ). The 
leading target mass corrections are given by the second term. As is discussed 
119, 120, 121, 122], these kinematical power corrections oc (M 2 /Q 2 ) n 



m 



may be summed up to all orders by the use of Nachtmann moments ||104j] . 
The third moment of G\ in the mass term of Eqn. (|54|) can be determined 
from data. The integral amounts to 0.0168 at Q 2 = 2 GeV 2 and 0.0130 
at Q 2 = 4.6 GeV 2 ||123|| . Target mass effects are therefore negligible for 
experimentally accessible Q 2 . One thus arrives at [35 



5S 



Bj 



(0.015 ±0.013)GeV 2 

Q 2 



(55) 



This result is comparable with a diquark model estimation of [124]. 

Higher twist effects were also studied for neutrino-nucleon scattering. In 
[ |125|| both the Baryon sum rule (see also [ |126|| ) and the Second Isospin 
SUM rule are discussed. With the reduced matrix element 



(P\0% m \P) apiD =2P,((Ol /m )) 

averaged 

of the twist four operators operators 



(56) 



(57) 
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one has 



5S t 



GLS 



5S™ poL 



3 ~ 



((Op) 

21) g 2 

((^5 NS )) 



(58) 



97 Q 2 

Target mass corrections may again be summed to all orders through the use 



of Nachtmann moments as discussed in [125] and are very small. The values 
for the reduced matrix elements were estimated as | |125 



{{Of)) = 0.33 GeV 2 
((Of 9 )) =0.15 GeV 2 



(59) 



Within the large uncertainty of ~ 50% this value is in agreement with the 
results from two other methods also studied in ||125|| , namely the vector dom- 
inance approximation and the nonrelativistic quark model. This indicates 
problems in bag-model calculations ||127|| which give negligibly small values 
for twist four contributions. The reanalysis of ||117|| presents as a result 
((Of)) = 0.53 GeV 2 with an error of about 20%. 

We are now in a position to compute the numerical values for the sum 
rules and their various corrections. As input parameter we take A 



(5)_ 
MS 
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MeV, where the superscript denotes the number of active flavors. This value 
corresponds to «( 5) (Mf) = 0.12. The running coupling constant a^J 1 (Q 2 , A^^ 
is calculated via the formula in the appendix. The matching equation at 

Q 2 



Ml 



a 



(M 2 A 



a 



^(M 2 ,A W 



1 



( n /) ' 

A/5 > 
(«/) 

a s 



^ H ' MS- 



TV 



(60) 



relates the coupling constant in the effective n/ — 1 flavour theory to the full 
theory with rif active flavors ||128| , |129| , |130|| . It is implicitly solved in or- 



Q 



M 2 



,(n f ) y 

(4.7 GeV) 2 and Q 2 



der to extract AjM . Applying this procedure subsequently at the thresholds 

MS 



Ml 



[1.6 GeV) 2 leads to Ai 



320 
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Q 2 


n f 




O(a s ) 


O(al) 




0{al) 


+ HT 






s G 


4 


4 


0.315 


0.333 


0.335 


- 


- 


_ 


5gls 


3 


3 


0.336 


3 


2.571 


2.583 


2.486 


2.388 


Sa 


3 


3 


0.336 


2 


2 


2 


2 


_ 


Qunpol. 
^Bj 


3 


3 


0.336 


1 


0.905 


0.902 


0.872 


0.828 




10 


4 


0.258 


1 


0.930 


0.931 


0.920 


0.907 




2 


3 


0.385 


0.210 


0.175 


0.175 


0.165 


0.157 




3 


3 


0.336 


0.210 


0.180 


0.180 


0.173 


0.168 




5 


4 


0.299 


0.210 


0.184 


0.185 


0.181 


0.178 




10.4 


4 


0.256 


0.210 


0.188 


0.189 


0.186 


0.185 


qp 
D EJ 


3 


3 


0.336 


0.185 


0.165 


0.166 


_ 


_ 




10.4 


4 


0.256 


0.185 


0.172 


0.174 






Cn 
°EJ 


3 


3 


0.336 


-0.024 


-0.0147 


-0.0146 








10.4 


4 


0.256 


-0.024 


-0.0154 


-0.0154 






<SbC 


3 


3 


0.336 


















10 


4 


0.258 

















Table 1: Theoretical estimates of the various sum rules. Corrections are 
taken into account up to and including the order indicated in the columns. 
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MeV and A^% = 366 MeV respectively. The various corrections to the differ- 
ent sum rules are displayed in Table 1. For the contributions of orders a s , a s 
and a z s we used the the formula for the effective coupling constant including 
the leading, next-to-leading and next-next-to-leading terms respectively. 



3 Comparison of Theory and Experiment 

3.1 Experimental Issues and Results 

Experiments have been performed at the Stanford Linear Accelerator Cen- 
ter (SLAC) using polarized and unpolarized electron beams, at FermiLab 
near Chicago using neutrino and unpolarized muon beams, at the European 
Center for Nuclear Research (CERN) in Geneva using polarized and unpo- 
larized muon beams and neutrino beams and at the Deutsches Elektronen- 
Synchrotron (DESY) in Hamburg using electron beams. Experiments are 
performed over a restricted range of x and Q 2 . Since the QCD corrections 
to the sum rules depend on Q 2 , data are required over the complete range 
of x in as narrow a Q 2 range is practicable. The range of x is restricted 
to x > Q 2 /s, In order for the QCD Parton Model to make a reliable pre- 
diction Q 2 ^,2 (GeV) 2 , hence, for the sum rules to be measured, data must 
be extrapolated into the very small x region. The extrapolation is least in 
experiments at the highest energy. To illustrate the extrapolation, consider 
the baryon SUM rule measured in neutrino scattering. Figure 2 shows the 
CCFR data ||131|. As well as F 3 , it shows / . dxF s (x) as a function of x r 



j mim 



the lowest value of x where data are available is x = 0.015. In performing the 
extrapolation to x m i n = a form for xF^{x) must be assumed. A fit of the 
form xF 3 (x) = Ax b (l — x) c provides an excellent description of the data. A 
systematic error must be included in the quoted value of the sum rule to take 
into account the extrapolation to x m i n = 0. This systematic error is difficult 
to estimate, since there is no fundamental reason for preferring one extrapo- 
lation over another. Data from different values of Q 2 can only be combined 
if a Q 2 extrapolation is assumed. Such an extrapolation can be based on a 
fit to perturbative QCD. However, if this is done, a "test of QCD" from the 
sum rules is compromised since perturbative QCD necessarily restricts the 
values that the sum rules can take. 



25 



GLS gum Rule: CCFR Data at Q 2 = 3 GeV 2 
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Figure 3: The Baryon sum rule from the CCFR experiment as a function 
of Q 2 . The inner error bar shows statistical errors and the outer a combina- 
tion of statistical and systematic errors. The curves are a QCD predictions 
that are discussed in the text. 



Having given these caveats, we will now discuss the current experimental 
values for the sum rules. For the baryon sum rule, the CCFR collaboration 
EH gives, 

S GLS = 2.50 ± 0.018(s*at) ± 0.078(sys) (61) 

at Q 2 = 3 GeV 2 . The systematic error includes the error from extrapolation 
into x = 0. A more precise result has been obtained by the CCFR col- 
laboration | 132|| by combining their data with that from other experiments 
on neutrino scattering ||133j , [134| , |135| , |136| , |137|| . Also, in the very large x 



region the nucleon's antiquark content is negligible and a relation between 
and Fl v exists in the Parton Model, viz 5(F 3 ^(x, Q 2 ) + F» n (x,Q 2 )) = 
18(F 2 ep (x, Q 2 )+Fg n (x, Q 2 )). The experimental precision on the baryon SUM 
RULE accuracy can therefore be improved by including data from F^ p ^ [ |138 



which has much higher statistics. The combined data are then extrapolated 
into the region below x = 0.02 and the sum rule evaluated. Figure 3 shows 
the extracted value of the baryon sum rule as a function of Q 2 . The 
curves on this figure will be discussed below. The solid line shows the QCD 
prediction including the higher twist effects, the dashed line shows the pre- 
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diction of the higher twist terms are ignored. A comparison of the two curves 
shows that the higher twist contributions are unimportant for Q 2 ^5 GeV 2 . 
The QCD prediction, which corresponds to a s (M z ) = 0.12, lies considerably 
below the data. The CCFR collaboration has fitted to QCD by allowing a s 
to vary ||132| (see also ||126| , [139| , |140|| ) . The best fit corresponds to 



a s {M z ) = 0.108 ± 0.004(stat) ± 0.004(syst) ± 0.006(FT) (62) 

The data for F% p — F% n and the integral f x (F2 P — Fffidx/x are shown in 
Figure 4 from ||141|| . The experiment observes no significant Q 2 variation over 
the range 0.5 GeV 2 < Q 2 < 10 GeV 2 . The measured values are extrapolated 
into x = assuming F 2 ep — F| r 



= ax . a and b are determined from a fit 
in the region 0.004 < x < 0.015 to be a = 0.2 ± 0.03, b = 0.50 ± 0.06. 
This extrapolation then contributes to the quoted error. The value of the 
Valence Isospin sum rule is determined by the NMC [ |141| | collaboration 
to be 

S G = 0.235 ± 0.026 (63) 

for < Q 2 >= 4 GeV 2 . This value is shown on the figure. The same experi- 
ment has issued preliminary results from its full data set [ |142| | which extends 
to smaller values of x and has fitted its data together with that of BCDMS 
143|| and SLAC ||138|| to give values of F^ and F| d (d represents deuterium) 



over the x range 0.006 < x < 0.9. Q corrections are applied to take higher 



twist effects into account and the results can then be interpreted as [144 



S G = 0.216 ±0.027 



(64) 



for < Q 2 >= 4 GeV 2 with no significant remaining Q 2 dependence in the 
range 0.5 < Q 2 < 10 GeV 2 . Recent data from E665 ||| agree with NMC. 
The Call AN- Gross relation is poorly determined. Experiments measure 



R(x,Q 2 



F 2 (x,Q 2 )(l + AM 2 x 2 /Q 2 )-2xF 1 (x,Q 
2xF 1 (x,Q 2 ) 



(65) 



Data from SLAC ||146| , |147| , |148|| are shown in Figure 5. The figure shows that 
R(x, Q 2 ) falls rapidly as Q 2 is increased and that it is small. The value is 
consistent with that predicted by QCD. The Isospin and Second Isospin 
SUM rules, which are related by the Callan- Gross relation, are difficult 
to measure with precision as they require neutrino scattering off hydrogen 
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0.005 0.010 



0.050 0.100 



0.500 1.000 



Figure 4: Data from the NMC collaboration on the VALENCE ISOSPIN 
SUM RULE. Shown is the quantity F| p — F| n (open points, right scale) and 
fx(F% p — F% n )dx/x (crosses, left scale). The bar at the extreme left shows 
the derived value of the sum rule. 




Q a (GeVY 



Figure 5: The quantity R(x, Q 2 ) from the SLAC experiments as a function 
of Q 2 for certain values of x at x = 0.5 (squares), x = 0.35 (diamonds) and 
x = 0.2 (crosses). 
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and deuterium targets and the statistical errors on such measurements are 
poor. Data show no significant Q 2 variation in the range 2 GeV 2 ^Q 2 ^30 
GeV 2 and give 

S A = 2.02 ± 0.40 (66) 

No published results are available on the Second Isospin sum rule. How- 
ever the result of Eqn. (|6~6"D together with the information on R(x, Q 2 ) show 
that the value of this rule is consistent with the expectation of S^" 1301 ' = 1. 

In the case of the Spin sum rules 
CERN " 



150, 151 



data are available from SMC at 
using polarized muon beams scattering off deuterium and 



hydrogen targets, and E142 ||152|| and E143 [153, |154|| at SLAC using polar- 
ized electron beams on He 3 , hydrogen and deuterium targets targets. The 
SLAC data only cover x > 0.03 and 1 GeV 2 £Q 2 £6.5 GeV 2 , while the CERN 
data extend to x = 0.003 and 1 GeV 2 ^Q 2 < ^60 GeV 2 . The experiments 
actually measure the asymmetry in scattering i. e. the ratio of the quantity a 
of Eqn. (|) to the unpolarized rate of Eqn. ([[]). No Q 2 dependence is observed 
in this ratio. 

In order to extract Gi(x, Q 2 ) it is assumed that this ratio is independent 
of Q 2 and therefore that the Q 2 dependence of G±(x, Q 2 ) is given by that 
of Fi{x,Q 2 ) and F 2 (x,Q 2 ). Gi(x,Q 2 ) is then extrapolated to Q 2 = 5 GeV 2 . 
Data are extrapolated to x = using the assumption that Gi(x, Q 2 ) ~ x a 
with < a < 0.5. An extrapolation in the region 0.6 < x < 1 is also needed, 
but this introduces a very small error since Gi(x,Q 2 ) is very small in this 
region. The contribution of the structure function G 2 (x, Q 2 ) is suppressed 
by the 1/Q 2 term in Eqn.@ and no information about it can be extracted 
from the data. 

These extrapolations enable the Spin sum rules for Gi(x,Q 2 ) to be 
evaluated. The SMC data alone ||150| , |I5T| , |155|| give for the Polarized 
Isospin sum rule 

S Bj = 0.199 ±0.038 (67) 

at Q 2 = 10.4 GeV 2 , whereas at Q 2 =3 GeV 2 , the SLAC data JT52|, |T53|, [T5i| 
give 

S B j = 0.149 ±0.014 (68) 

The different form used for the extrapolation to smaller values of x is partly 
responsible for the smaller values. The experiments can be combined with 
earlier results involving hydrogen targets ||156j, |157|| to give, at Q 2 = 5 GeV 2 , 
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o 2 

OeV 2 


theory 


Expt. 


s G 


4 


0.335 


0.216±0.027 


Sgls 


3 


2.388 


2.50 ± 0.08 


Sa 


3 


2 


2.02± 0.40 


Sbj 


5 


0.178 


0.203± 0.023 


CP 
°EJ 


5 


0.171 


0.136± 0.010 


On 


5 


-0.0135 


-0.067 ± 0.016 



Table 2: Comparison of the theoretical and experimental values of the sum 
rules. 



S£j= 0.136 ±0.010 

(69) 

S£ J= -0.067 ±0.016. 
The Polarized Isospin sum rule determined from these is 

5 Bj = 0.203 ± 0.023 (70) 

3.2 Theory vs. Experiment 

Table 2 shows a comparison of the experimental values discussed above with 
theoretical predictions from Table 1. In the cases where the experiments have 
corrected for the effects of higher twist the relevant comparison is with the 
highest order perturbative QCD result available and it is this number that is 
given in the theory column. No entries are shown for the G2 and Second 
Isospin sum rules where no data exist. 

It can be seen from the Table that the sum rules fall into three categories. 
First, the Isospin sum rule has very large experimental uncertainties but 
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the measured values are consistent with the expectations of QCD. Second, 
the Baryon sum rule and Polarized Isospin sum rule are compatible 
with QCD, but have experimental errors that are small enough so that the 
measurements can discriminate between the QCD results at different orders 
in perturbation theory. In these cases the data are consistent with the QCD 
expectations and are inconsistent with the Naive Parton Model. Finally, the 
Spin and Valence Isospin sum rules have experimental values that are 
inconsistent with the Naive Parton Model or QCD predictions. 

The second category can be used to measure the strong coupling constant 
a s . Figure 3 shows the Q 2 dependence of the Baryon sum rule. This 
value is somewhat lower than the world average [|159|| The quoted error is 
dominated by that due to the Higher Twist terms (HT). The Polarized 
Isospin sum rule has also been used to determine a s [|160|| :- 



a s (M z ) = 0.1221^09 (71) 
if the higher twist terms are neglected, or 

a s (M z ) = 0.118+°;°?! (72) 

if they are included. 

The final category needs more discussion. The Valence Isospin sum 
RULE discrepancy between theory and experiment shown in table 2 can be 
removed by dropping the assumption that J 1 dx(u(x) — d(x)) = 0. Using the 



Naive Parton Model (see Eqn. [T4|) and the data we obtain 



f dx(u(x) - d{x)) = -0.176 ± 0.040. (73) 

J 

The QCD corrections are much smaller than the error on this result. Addi- 
tional experimental information is available from the processes pd — > n + fi~ + 
anything and pp — > /i + yU~ + anything which, in the Parton Model, are 



due to quark antiquark annihilation. Data from NA51 ||161|| indicate that 
u(x)/d(x) = 0.51 ± 0.08 at x = 0.18 and Q 2 = 2.5 GeV 2 . The possible 
non-equality of u(x) and d(x) was first suggested in ||162|| where a possible 



parameterization was introduced. Several authors have attempted to esti- 
mate the size of u(x) — d(x) that could arise from non-perturbative effects. 
Some have attempted to explain the effect in terms of a pion cloud surround- 
ing the nucleon [ |163|| . Other models are based a chiral Lagrangian ||164 
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approach that starts with a nucleoli consisting of three valence quarks then 
generates the anti-quark distributions from pions emitted in processes like 
u — > d-n + | |165| , |166| |. Models of this type produce a difference (u(x) — d(x)) 
that is concentrated at very small values of x and the ratio u(x)/d(x) is 
predicted to be quite small. 

The predictions for the values of the Spin sum rules depend upon the 
assumed values for ao, a 3 and a%. There is no fundamental reason for the first 
of these to take the value ao = v^3a§ which was the original assumption of 
[ 23|| and is the value used in the table. If we use ao = 0.12 ± 0.027, obtained 



from the value of AS" obtained from up elastic scattering ][|] we have at 
Q 2 = 5 GeV 2 the "predictions" Sfa = -0.065 ± 0.030 and S& 3 = 0.12 ± 0.03 
both of which agree with the experimental results shown in Table 2. Instead 
we can use the experimental results for S P EJ and S EJ to determine ao from the 
QCD forms of Eqn.(|4*2"D together with the higher twist corrections. This gives 
ao = 0.19 ± 0.07. Using the Naive Parton Model relation of Eqn.fllOD implies 
that J dx(As(x) + As(x)) = —0.13 ±0.03 from whence we can also infer that 
J dx(Au(x) + Au(x)) = 0.8 and / dx{Ad(x) + Ad(x)) = 0.4. Everything is 
consistent but one is left with the annoying question of what is carrying most 
of the nucleon's spin. In the model where the nucleon is viewed a soliton-like 
solution of ||167|| one expects ao = ||168j| in the limit of zero quark mass 



and large number of colors. In this interpretation all of the nucleon's spin is 



carried by orbital angular momentum. Reference [|169|| can be consulted for 
a detailed review. 

If the gluon distribution in the proton is polarized, there is an additional 
complication. The scattering process e + g — > eqq can generate a contribution 
to Gi(x, Q 2 ) at order a s . Adding this term to the Naive Parton Model result 



is equivalent to replacing Eqn.(|l0|) by ||170 | 



a 



J dx(Au + Au + Ad + Ad + As + As - nf^Ag^j (74) 



This contribution is not present in the Operator Product analysis presented 
above. It can be introduced if one observes that while the operator cor- 
responding to the singlet axial current is not conserved and is therefore 
subject to renormalization due to the axial anomaly [ 171 , |172j , 173 ], a lin- 



ear combination of this operator and a gauge variant operator made up of 
gluon fields is not renormalized. If this term is included the form of the 
QCD corrections given in Eqn.(f4"2"|) are the same except that a nv is inter- 
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preted as that of Eqn.([74]). The data are now to be interpreted as implying 
fdx(As + As - n f ^Ag = -0.13 ± 0.03. If we assume that AS = then 
/ dxAg = 3.2 ± 0.75 at Q 2 = 5 GeV 2 . This substantial polarization should 
be observable in other experiments. For example, the production of pions at 
large transverse momentum in proton-proton scattering proceeds via parton 
parton scattering of the type q + g — > gq. If both protons are polarized an 
asymmetry 



A 



LL 



- da(-+) 

da(++) + da(-+) 



(75) 



can be formed (the ± arguments refer the helicity of the incident protons) 
which is depends upon Ag(z). An experiment at FermiLab [|174j| observes 
an asymmetry that is consistent with zero for transverse momenta of pions 
less than 3 GeV. More recently ||175|| the same experiment has measured the 
asymmetry for double 7r° production. Again the asymmetry is consistent 
with zero. If models for Aq are assumed ||1 76|| then a constraint can be 



obtained on Ag. This constraint is sufficient to rule out some models ||177 
but others that have / dxAg(x) ~ 5 JT7BJ are not excluded. 



4 Conclusions 

The Parton Model sum rules represent fundamental predictions of QCD. The 
experimental precision of many of these rules is such that consistency with 
the theory can be established. In the case of a few of the rules, notably 
the baryon SUM rule, the data are sufficiently precise that consistency 
can be checked in detail and a value of the strong coupling constant obtained 



whose error is competitive with the best measurements [|159|| . In this case the 
theoretical errors coming from the poor knowledge of higher twist terms and 
the order terms contribute significantly to the error on a s . Improvement 
in these areas is unlikely to appear in the near future. The failure of the 
Valence Isospin sum rule has led to the realization that u(x) ^ d(x) 
and while there is some theoretical understanding of how this might arise, 
the difference, like all other structure functions, must be extracted from 
data. The failure of the naive form of the Spin sum rules has led to an 
interesting situation. There must be significant polarization in the strange 
quarks and/or the gluons. More accurate data on up elastic scattering might 
enable the former to be constrained. The latter should be constrained when 
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polarized proton-proton scattering experiments become available at RICH in 



the next few years [179 



The advent of data from HERA JTHOfl have enabled structure functions to 
be measured at smaller values of x than those in fixed target experiments. 
Nevertheless, the statistical errors on these data are still quite large and they 
are, of course, only available for F% p (x, Q 2 ). In the future, data from polarized 
ep scattering will be available from this facility ||181 that will considerably 
extend the range of x and Q 2 available for the measurements of G\(x, Q 2 ) and 
reduce the error on the Spin sum rules resulting from the extrapolation 
into x = 0. 
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Appendix 



The effective coupling constant in next-next-to-leading order may be written 
in the form 



71 



1 



1 PilnL 



m U 2 



(In 2 L - In L - 1) + 



02 



(76) 



where L = ln(Q 2 /A 2 ) and a s {Q 2 = fi 2 ) = a s (fi 2 ). The coefficients of the 
beta-function are known up to the three loop level 0, ^, 182 , 183 , 184 , 185 
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Anomalous dimensions of singlet and nonsinglet operators were calculated 
in a number of works for both unpolarized [ |59l , 60, 61, 32], [39, |63|, 34, 
|6T| [67| and polarized scattering [52], [53], [>|, [5^, Some 



of the results are given here. The coefficients of the nonsinglet anomalous 
dimension for unpolarized scattering read 



NS(0) 
m=l 



NS(1)- 
ln=l 



NS(1)- 
in=l 



4 



n(n + 



n=l 



--(13 + 8C(3)-2tt 2 
0. 



(78) 



For polarized scattering one has the following nonsinglet anomalous dimen- 
sions A 7 ^ s (°) = 7^ s (°) and A^ 1 ^ = 7^ s(1)± . Finally the coefficients for 
the singlet quark diagonal anomalous dimension matrix elements are 



A7^ (0) 



A 7; 



S(l) 
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A 7 



S(2) 
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Finally we give the analytic formulae of the coefficient functions corre- 
sponding to the first moments of the various structure functions. The nonsin- 
glet coefficient function for the structure function F 2 of unpolarized electron 
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nucleoli scattering is given in the MS scheme |)9|, [62|] by 
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(80) 



The nonsinglet JT^ |0|, 0, ||, and the singlet @, |3|, |3| coefficient 
functions for polarized electron-nucleon scattering read 
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The coefficient function for the neutrino structure function Fa has the fol- 



lowing form H, H, H, ||]: 
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The coefficient function for the nonsinglet structure function F\ of neutrino 
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and antineutrino scattering reads p3, p3, 190 
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